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SUMMARY 

An approximate algebraic technique which extends the usefulness of three-gimbal 
systems is presented. In analyzing the three-gimbal problem, the Euler angles are 
explored about the singular points and are found to be, in general, non-unique. This 
non-uniqueness is eliminated by minimizing the magnitude of the discontinuities at the 
singular points. An approximate mathematical model of the inverse Euler angle prob- 
lem is developed which chooses the solution with minimum discontinuity and, in addition, 
replaces high rates over large angles by small discontinuities. A computer circuit is 
devised for this model. 

A proposed four-gimbal system is investigated and a mathematical model is derived 
by finding a necessary differential constraint on the fourth angle to maintain a nonsingular 
formulation. This constraint is then used to derive a particular mathematical formulation 
in an essentially algebraic form, from which a computer circuit is developed. The rate 
equations of the four angles are also derived. Existing four-gimbal stable platforms are 
analyzed and discussed. 


INTRODUCTION 

Three-gimbal systems have been used for several years as a means of displaying 
visual information in flight simulation. Traditionally the Euler angles, needed to drive 
these gimbals, were also used in the main calculations as parameters of rotation from 
which the direction cosine matrix was computed. Difficulties with the singular points 
were normally avoided since only small portions of the flight program were studied. 
Thus, gimbal systems could be chosen with gimbal lock out of the region of interest. 

As the need to study highly maneuverable fighter aircraft and spacecraft arose, the 
trend was to remove the errors due to singularities from the main calculations by using 
either the Euler parameters or direction cosine rate equations. The Euler angles were 
then determined algebraically from the direction cosines (ref. 1, pp. 415-418). Although 



the main calculations did not suffer from the singularities, little concern was given to the 
resultant positioning errors of the visual display system (ref. 2). 

The newest concept in flight simulation is the full-mission simulator. This con- 
cept has had a tremendous impact on both simulation computing practices and visual dis- 
play hardware design. Present-day full-mission simulators require complete freedom 
of rotation of the coordinate axes with the ability to hold at any altitude for indefinite 
periods and angular drives for several orders of rotation between several different such 
coordinate axes. To meet these requirements, an algebraic computation is a must in 
order to synchronize the motion, and considerable care must be maintained in the com- 
putation of the transformation matrices between these various axis systems. Still an 
open question is the method to determine the Euler angles in the neighborhood of the 
singular points. 

The behavior at the singular points is not a serious consequence until mechanical 
equipment (namely, a three gimbal) is required to track the motion. Then consideration 
should be given to minimizing the dynamic positioning error due to rate limiting of the 
gimbal system. The purpose of this report is to develop a simple method by which the 
positioning errors can be reduced. In minimizing the dynamic positioning error, it is 
found that a static error is introduced in the neighborhood of the singular points. This 
method contains an adjustable parameter in the form of a dead band about the singular 
points. A value for the dead band, within computer capabilities, can be chosen which 
best meets the mission requirements without driving the gimbal far beyond its response 
limitations . 

Even though the method established in this report can reduce the positioning errors 
associated with the singularities of the three-gimbal system, some rate limiting and static 
error still exist in the neighborhood of the singular points. A four -gimbal system can 
remove these difficulties. In this paper are presented the development of a mathematical 
model for a specific four-gimbal system and the results of simple experiments to demon- 
strate its operation. A reasonable amount of ingenuity should allow extension of this 
analysis to any general four-gimbal configuration. 

With the increased capability of modern aircraft and spacecraft, flight situation 
indicators that can adequately display attitudes of the vehicles are required. To meet 
this need, four-gimbal stabilized platforms have been proposed (refs. 3 and 4). How- 
ever, the methods proposed to drive these systems yield discontinuities similar to those 
encountered in three-gimbal systems. When these discontinuities occur, the corre- 
sponding gimbal drives must be delicately balanced to assure platform isolation. These 
systems are analyzed herein and the cause of these discontinuities is ascertained. 
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SYMBOLS 


Function: 

A(t) transformation matrix from fixed- coordinate frame to rotating- coordinate 

frame, dimensionless 

ay elements of A(t), where subscripts i and j = 1, 2, or 3, dimensionless 

a,b,c,d Euler parameters of rotation of rotating-coordinate frame with respect to 

fixed-coordinate frame, dimensionless 

C y cosine of any angle y, dimensionless 

K positive gain constant 

k,Z,m,n integers 

p,q,r components of angular velocity vector of rotating-coordinate frame with 

respect to fixed- coordinate frame, measured along rotating coordinate 
axes, radians/second 

Sy sine of any angle y, dimensionless 

sgn( ) signum function, equal to +1 for variable 50 and -1 for variable <0, 
dimensionless 

t time, seconds 


Vf vector measured in coordinate frame considered fixed, arbitrary units 

V r vector measured in coordinate frame considered rotating, arbitrary units 



Cy Sy 0 
-Sy Cy 0 
0 0 1 
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tr\ 


inverse or transpose of fyj , dimensionless 

A 


Euler angles relating fixed frame to rotating frame in three- gimbal system 
(bank angle, pitch angle, and heading), radians 

a ,$,#,<// angles relating fixed frame to rotating frame in four-gimbal system (for a = 0, 
<(> is bank angle, 9 is pitch angle, and <// is heading), radians 

maximum slew rate of gimbal system, radians/second 

cu angular velocity vector of rotating- coordinate frame with respect to fixed- 

coordinate frame, radians/second 

a T q a; measured in three-gimbal space, radians/second 

Subscripts : 

o value at t = 0 

+ signifies right-hand limit 

signifies left-hand limit 

Dot over symbol indicates derivative with respect to time. 

Computer: 

servomotor that turns shaft in direction of input voltage 

integrating servomotor that turns shaft proportionally to integral of 
input voltage 





potentiometer whose output voltage is K times input voltage 
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inverting amplifier whose output voltage is the negative of input voltage 


e i j 

ABS 

absolute value circuit whose output voltage is equal to magnitude of 
input voltage 

iV, 

op. 

1 + 

relay such that, if the sum of input voltages is negative, the switch is in 


T- 

the - position and in the + position otherwise 

i e i 

l e o 

y<— 

shaft driven resolver whose output voltages are 

l e o = cos 9 - 2 e i sin 0 

2 e i 

M 2 e o 

2 e 0 = 2 e i cos e + l e i sin 0 


f. 

where 9 is shaft rotation 



THE THREE-GIMBAL PROBLEM 



Analysis 

The order of the Euler angles (fig. 1) considered in this analysis is that of the con- 
ventional airplane angles if/, 9, and $. It is assumed that the transformation matrix 
A(t) is available as a function of time and is defined by 


V r = A(t)V f 

The transformation matrix was generated by the Euler param- 
eters given by equations (Al) and (A2) in appendix A. 

The transformation matrix A(t) uniquely defines the 
relative orientation of the two coordinate frames as a function 
of time. The Euler angles are then related to the elements of 
A(t) as 

s ^ ( 2 ) 

The inverse Euler angle problem is to solve these functional 
relations for the three angles (that is, mapping of the elements 
of A(t) into ip, 9, and rp). Equation (2) is expanded in 


( 1 ) 



Figure 1.- Three-gimbal configuration 
considered in the analysis. 



equation (A3) of appendix A. The inverse problem is seen to be everywhere non-unique * 
since 

1$ ± 4J? ~ ± = 


The inverse of the mapping exists and is locally one-to-one when the Jacobian 
does not vanish (ref. 5, pp. 225-232), otherwise it is said to be singular. As cos 9 
approaches zero, the Jacobian of the inverse problem vanishes; that is, no three ele- 
ments of A(t) can be found such that 


9 ( a ij> a kl> a mn) | 

a 0 P,o,4>) 


* 0 


cos0=O 


where 

i,j,k,l,m,n = 1, 2, or 3 

Hence, a functional dependence exists between if/ and cj>. This functional relation is 
found from equation (2) as cos 9 approaches zero and sin 9 approaches ±1 to be 


0 0 Tl 

±S(<f> T if/) C(cp T if/) 0 
±C(<p T if/) -S((p T if/) 0 


= A(t) 


(4) 


An alternate method of computing the Euler angles is the following Euler angle rate 
equations (ref. 1, p. 351): 

rC(p + qS <f> 


* = 


C9 


(5) 


9 = q C<p - rSc p 


( 6 ) 


4> = p + ij/S9 (7) 

The solution to the rate equations is unique on intervals not containing the singular 
points. If a singular point is encountered, uniqueness is no longer assured. It is shown 
later that, as cos 9 approaches zero, equation (5) is indeterminate. At the singular 
points, sin 9 = ±1, equation (7) can be rewritten (ref. 6) as 

<p t P - j* P dt (8) 


^Non-uniqueness is restricted herein to solutions within ±tt of the principal 
solution. 
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Indicative of both the inverse problem and the rate formulation at the singular 
points is the functional dependence of if/ and cp; these angles are not defined at the 
singular points but either the sum (at sin 9 = -1) or difference (at sin 9 = +1) is 
defined. On every interval not containing a singular point, each solution of the inverse 
problem satisfies equations (5) to (7). The solution to the rate equations cannot be 
propagated uniquely through the singular points, but every set of functions which satisfy 
equations (5) to (7) except at the singular points will correspond to the solutions to the 
inverse problem. Hence, the analysis of one formulation discloses the behavior of the 
other. 


The remainder of the analysis relies heavily on the principle of invariance. A 
transformation (of a vector) of invariance is such that if the components of the vector 
are transformed invariantly, then the vector is unchanged (that is, the transformation 
does not alter its magnitude or direction). 

Consider the space with three coordinate axes defined by the directions of increasing 
< P , 0, and shown in figure 2. These directions are the same as those of the gimbal 
axes shown in figure 1. The angular velocity vector uq, referenced to this oblique 
three-dimensional space, is given by 

w G = TcJ (9) 


The angular velocity vector co is referenced to a three-dimensional orthogonal space 
with components p, q, and r. The transformation T transforms from the orthogonal 
three-dimensional space to an oblique three-dimensional space with coordinates \p, 9, 
and as shown in figure 2. Using equations (5) to (7) with equation (9) gives 


0 

Sep 

Cep 

ce 

C9 

0 

Cep 

-Sep 

1 

se Sep 

S9Cep 

C9 

ce 


( 10 ) 


As seen from figure 2, the obliqueness of the space 
depends on the angle 9. In the limit as 1 9 | — tj/2, the 
oblique three-dimensional space collapses into a two- 
dimensional orthogonal space which is capable of repre- 
senting only angular velocities that lie in the plane of 
the space; only the projection of a general vector in the 
three-dimensional space onto the two-dimensional space 
can be represented as seen in figure 3. (Note that the 
two-dimensional space coincides with the plane of the 
gimbal at gimbal lock.) 



Figure 2,- Ordinary three-dimensional space and 
oblique three-dimensional space related through 
transformation T. 
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If invariance of the angular velocity vector under change 
of representation is required for every value of 6 except for 
1 0 1 = ir/2, then the angular velocity vector must lie in the plane 
of the space (that is, in the limit as | 0 1 approaches v/2). In 
the limit as the singular point is. approached, the out-of-plane 
component must vanish. The left-hand limit (time increasing) 
and the right-hand limit (time decreasing) are, in general, not 
equal (that is, the solution is discontinuous). This limiting pro- 
cess can be loosely viewed geometrically; if an angular velocity 
vector "appears" that does not lie in the planar space, then the 
space must rotate instantaneously to contain the vector. 

The qualitative behavior at the singular points is established. In order to look at 
the limiting process, the problem is reformulated explicitly in terms of the out-of-plane 
and in-plane components by rewriting equations (5), (6), and (7) as 


$ - }J/S6 


O 

o 
1 ^ 

p 

9 

= 

0 C (p -S <p 

q 



0 S (f> C(p 

r 


The vector on the left-hand side of equation (11) is referenced to a three-dimensional 
orthogonal space. Let the components be x', y’, and z', where the z’- component is 
recognized as the out-of-plane component in the limit as indicated in figure 4. 



Figure 3.- Oblique three-dimensional 
space collapsed to a two-dimensional 
orthogonal space at singular points. 
The in-plane and out-of-plane com- 
ponents of a general vector in three- 
dimensional space are shown. 


Figure 4 shows the three-dimensional space with 
components x', y', and z’ (at the singular points) in 
relation to the three-dimensional space with components 
p, q, and r to which To is referenced. The transfor- 
mation from p,q,r to x'jy^z’ is given by 


x T 


O 

O 

p 

y* 


0 Cep -&<p 

q 

z' 


0 S <f> Cep 

r 


( 12 ) 


±<l> 

X' 



At the singular points, the X’Y' -plane is the planar 
space of the gimbal shown in figure 3. The left-hand limit 
of the solution <p as the singular point is approached is 
shown as the angle between the (Y’,0) and the q axes. 
After the singular point is encountered there are two solu- 
tions corresponding to (y',0)j and (y* , 0)2 for which 
the right-hand limit (time decreasing) of the out-of-plane 



Figure 4.- The two-dimensional space 
(XT-plane) is shown in relation to 
pqr space. The z'-component is the 
out-of-plane component that must be 
reduced to zero by the indicated dis- 
continuous rotations. 
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component vanishes, where of in figure 4 is the angular velocity vector at the exit (right- 
hand limit) from the neighborhood of the singular point. The right-hand limit of the out- 
of-plane component of the angular velocity vector is 

z' = q + sin cp + + r + cos <p + = 0 (13) 

which is the numerator of equation (5). Hence, \p is indeterminate. Equation (13) has 
two principal solutions corresponding to A and A </>2 in figure 4, which are the 
angles the space must be rotated instantaneously to contain the angular velocity vector. 
From figure 4, therefore, 

A0j = u - A$2 (14) 

As shown by equation (14), one discontinuity always has a magnitude less than or equal to 
w/2. As indicated by equations (4) and (8) a corresponding discontinuity occurs in \p. 

From this analysis, the following properties are noted: 

(1) The angles and cp are not defined at the singular points but either the sum 
or difference is defined. Hence, the angles can be chosen arbitrarily, subject to the pres- 
ervation of equations (4) and (8). 

(2) In general, discontinuities appear at the singular points but a discontinuity can 
be chosen to be n/2 or less in magnitude. 

Application 

A method is presented for computing the Euler angles from the elements of A(t). 
The Euler angles are functionally related to A(t) by equation (2). The non-uniqueness 
shown by equation (3) is used in minimizing the discontinuities at the singular points. 

The angle (p can be computed from a 2 g and agg (eq. (A3), appendix A) by using a 
bootstrap (iterative) approach which makes cos 6 available as 

aoo 

sin 0 = ZZTe (cos 6 * 0) (15) 

a 33 

cos 0 = ^Te (cos 9 * 0) (16) 

With this knowledge of <p, a new matrix B(t) can be defined as 

B(t) = [-<£] A(t) (17) 

A 
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This matrix is functionally related to ip and 6 by 

B(t) = (18) 

Equations (17) and (18) are expanded in appendix A. The solution of equation (18) is well 
defined, the form being that of the familiar two-gimbal problem. 

The difficulties are left in the solutions of equations (15) and (16) when cos 9 is 
near zero and in minimizing the discontinuities. To circumvent these difficulties, the 
two properties summarized in the analysis are used. 

A small dead band is established about the singular points within which 0 is held 
(by property (1)) at the value computed by equations (15) and (16) at the entrance to the 
dead region (the left-hand limit of 0). The computation of \p and 9 from equation (18) 
insures the preservation of equations (4) and (8) at the singular points. If a roll maneuver 
(p * 0) is initiated in the dead region, an effective roll occurs through the \ p angle and 
the visual display moves in the proper direction. 

As noted in the analysis, the discontinuities can be chosen to be tt/2 or less (prop- 
erty (2)). The magnitude of the discontinuities must be computed so that the proper choice 
can be made. To compute the discontinuities, the left- (time increasing) and right-hand 
(time decreasing) limits of the solution as the singular points are approached must be 
found. The left-hand limit 0 _ is computed by equations (15) and (16) upon entrance to 
the neighborhood of the singular points. With this value and the element b 33 of equa- 
tion (18), cos 9 is computed from 


^33 = cos 0 = a 23 s * n 0 + a 33 cos 4>_ 

(19) 

The current value of 0 is contained in a £3 and a 33 due to the following functional 
form: 

a 23 = cos 0 s * n < P + 

(20) 

a 33 = cos 9 cos 0 + 

(21) 

Hence, when 0 is discontinuous, equation (19) becomes 


bgg = cos 9 COS A 0 

(22) 

where 


A0 = 0 + - 0_ 

(23) 


is the discontinuity. 
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Maintaining 


I (24) 

provides the necessary choice between the solutions given by equations (2) and (3) (that 
is, cos 8 or cos(7r - 9) = -cos 9) which minimizes the discontinuities. Hence, the 
choice 

sgn(cos 8) = sgn(b 33 ) = sgn(cos 0)sgn(cos A $) (25) 

completes the necessary logic for a solution that minimizes discontinuities since equa- 
tion (25) implies that 

sgn(cos A <p) = +1 (26) 

which indeed reflects inequality (24). 

The dead band has the effect of introducing a static error and limiting the maximum 
angular rates of and 0 except, of course, for the discontinuities, the magnitudes of 
which will be minimized. Hence, there is a trade-off between static error and error due 
to rate limiting of mechanical equipment by choice of the size of the dead band. 

Mechanization 

Mechanizing the problem as formulated in the previous section can be accomplished 
in many ways. The method utilized herein uses the resolver chain described in refer- 
ence 7. This mechanization is shown in figure 5. In the figure, the dead band (DB) is 
detected by 

| a 23 1 + | a 33 1 = (l cos 4>\ + Jsin 0|)|cos 9\ 

which gives a measure of | cos 8\. Note that 
I a 13 I = | s * n ® | cou ld have been used as well. 

When in the dead region the 0 resolver is 
physically turned off by removing the power 
supply from the motor windings, as indi- 
cated in figure 5. The error signal to the 
< p resolver has its sign determined by 
sgn(cos 9 cos A$) which provides the mini- 
mization of discontinuities when the <p 
resolver is again made active. Inasmuch as 
the 9 and if; calculations are not directly 
affected by the dead band, equation (8) is 
preserved within a small error inside the 
dead region. 



Figure 5.- Euler angle computer circuit using resolver chain with asso- 
ciated logic for minimization of discontinuities. 
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A reasonable estimate of the size of the dead band can be made by observing the 
effects of the dead region on the maximum angular rates. The maximum angular rate 
of the system is determined from equation (5) which is 

r cos <b + q sin d> 
iL/ = 

cos 9 


Let the maximum slew rate attainable by the gimbal be Then it is desired that 

= W L 


from which (with eq. (5)) 


e db " 


COS 


sup(|q | + |r I) 


O). 


is the value for the dead band for which rate limiting will rarely occur except for the dis- 
continuities the magnitudes of which are minimized. This estimate does not consider the 
static errors inside the dead region. 


THE FOUR- GIMBAL PROBLEM 


Analysis 


The four-gimbal configuration is shown in figure 6. It is 
assumed, for the purposes of the following derivations, that 
the transformation of coordinates, matrix A(t), is available 
and defined as 

V r = A(t)V { (1) 

The transformation matrix A(t) can be constructed as a func- 
tion of the four angles \p, 9, <p, and a. This functional 

relation is 



A(t) = [a] z Efl x OOyl30 z 


(27) 


Figure 6.- Four-gimbal configuration 
considered in the analysis. 


and is expanded in equation (Bl) of appendix B. 

The introduction of the fourth angle makes it possible to move the singular points in 
the solution by the choice of the mathematical model. To demonstrate these singular 
properties, recall that singularities occur in the three-gimbal system when two axes are 
alined. Then several simple cases arise by alternately letting one of the four angles be 
constant. These cases are summarized in the following table: 
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Angle held constant 

Term causing singularity 

or = 0 

cos 0 = 0 

CM 

II 

a 

cos 0 = 0 

0 = 0 

sin 0 = 0 

-e- 

ii 

to 

Everywhere (two gimbal) 

0=0 

sin 0=0 

0 = 7r/2 

Everywhere (two gimbal) 

xp = 0 

cos 0=0 

xp = tt/2 

cos 0=0 


It is clear from viewing the table that many possible singularities can enter the solution. 
Formulations of the problem can be found that have singularities involving more than one 
angle, including the angles a and xp which do not appear in the right column of the 
table. The problem is simplified by retaining xp, 9, and 0 as a general rotation from 
a coordinate frame that is orientated by the angle a. The singularity then arises from 
the one term cos 0=0. The problem reduces to choosing a such that the singular 
point is never encountered. 

Assume that a value for oi has been chosen; thus far the choice is entirely arbi- 
trary since four parameters overspecify the rotation of two coordinate frames (ref. 8). 
With this knowledge of a, a matrix D(t) can be defined as 

°(t) = C-q] z A(t) = (28) 

These relations are expanded in equations (B2) and (B3) of appendix B. The solution for 
9, and 0 is singular at cos 0 = 0 and unique in the sense of principal values for 
cos 9 > 0. To remain an unique nonsingular solution for xp, 9, and <p, the angle a 
must be chosen such that cos 0 remains always greater than zero. 

The only difficulty left is the choice of a. Suppose that cos 9 is approaching 
zero; a must be changed by Aar such that the change in cos 9 (that is, A cos 0) is 
positive. The first-order change in cos 0 with respect to a is given by 

A cos 0 = - £° S 0 Act (29) 

and requiring cos 0 to increase yields 


Since 


9 COS 0 
9at 


Aar > 0 


Aar = ar At 


(30) 

(31) 


13 




Equations (31) and (32) result in 


L£OSJ_ 6l> o (32) 

da 

Once the partial derivative is evaluated, equation (32) can be used to determine the direc- 
tion to change a in order to increase cos 9. 

In order to evaluate the partial derivative in inequality (32), cos 9 must be 
expressed as a function of a, <p, and t. This functional form is given by the matrix 

C(t) = = [e] Y pO z (33) 

from which 

C33 = cos 6 = a^j sin <p sin a + a23 sin <p cos a + a33 cos <p (34) 

as seen in equations (B4) and (B5) of appendix B. Differentiating equation (34) explicitly 
with respect to a yields 

- C ^ QS — = a^'Sin <f> cos a - sin <p sin a = sin 0(ai3 cos a - a23 sin a) (35) 

From equations (B2) and (B3) of appendix B 

df3 = -sin 9 = aj3 cos a - a23 sin a 


and, hence, equation (35) becomes 


9 cos 6 = -sin 9 sin 4> 

da 

(36) 

From relation (32) 


(-sin 9 sin cp)a > 0 

(37) 


must be maintained to insure a nonsingular solution for i//, 9, <p, and a. Clearly 

a # 0 must be satisfied whenever cos 9 approaches zero. 

Relation (37) defines the direction a must be changed to increase (maximize) 
cos 9. Hence, the direction is defined by sgn(-sin 9 sin <p) and only the magnitude 
is left undetermined, that is, 

a = f(o!,i//,0,4>,t)sgn(-sin 9 sin 0) (38) 

where f(9) (in the sequel the arguments a, \}s, (p, and t are suppressed) is an arbi- 
trary nonnegative function which is nonzero for cos 9 approaching zero. A desirable 
form of f (9) can be chosen from two considerations: simplicity and the condition that 
a goes to zero smoothly as cos 9 approaches the maximum (that is, sin 9 approaches 
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zero). Such a function is 


f(0) = K|sin a| 

where K is a positive constant (to be determined) and this choice yields 

a = -K sin 9 sgn(sin <p) (39 

To investigate further the four-gimbal problem, the time derivatives of the four 
angles are needed. With \p f 9, 0, and a used as parameters of the rotation, the 
rate equations are derived by transforming the angular velocity components 9, <p, 
and a into the axes of the rotating frame. The resulting angular velocity vector com- 
ponents are then formally equated to p, q, and r; hence, 


q = 0 + [qQ z o + 9 + 0 


Expanding equation (40) gives 

fca i 


C$Sar 

C9S(j)Sa - S0Ca 

CffiCai 

C9S<pCa + S0Sa 

-S<f) 

C0C<£ 


P 

q 

r - a 


Solving equation (41) for if/, 6, and 0 results in 

(pSa + qCa)S0 + (r - a)C(p 


8 = (pSa + qCa)C(p - (r - a)S(j> (43) 

4 > = (pCof - qSo?) + S (44) 

which are of a familiar functional form and where a is given by equation (38) as follows : 

a = -f(0)sgn(sin 8 sin (f>) 

In the limiting case as a and a approach zero, equations (42) to (44) are the standard 
Euler angle rate equations- given by equations (5) to (7). 

Now consider the constraint inequality (eq. (37)) 

(-sin 6 sin <p)a > 0 
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Substituting equation (38) into equation (37) gives 


sin 9 sin 0|f(0) > 0 


( 45 ) 


If f(0) is to be bound for sin 6 * 0, then inequality (45) cannot be satisfied for 
sin 0 = 0. The problem arises whether or not sin 0 = 0 is a singularity. To deter- 
mine the behavior at sin 0 = 0, consideration must be given to the change in sin 0 due 
to variations in a. This change can be determined by explicit differentiation of d23 in 
equations (B2) and (B3) of appendix B where 


d £3 = cos 6 sin 0 = aj 3 sin a + a 23 cos a 

Differentiating gives 

. 9 sin 0 . , 9 cos 9 ■ n 

cos 9 — + sin 0 = a-n cos a - aoo sin a = -sin 9 

da da 

from which 

9 sin 0 -sin 6 <? , 

— = COS^0 

da cos 0 


The first-order change with respect to time of sin 0 due to & is 

37 (sin 0) = i(o) --- - - - -- cos20 sgn(sin 0) 

dt cos 9 


(46) 


which is nonzero when sin 9 and cos 0 are nonzero provided f(0) is also nonzero. 
The choice of f(0) to be nonnegative tends to maximize | sin 0 | as | sin 9 | approaches 
a minimum in such a way that the orientation of the coordinate systems is preserved. 
Hence, the choice of f(0) to be nonzero and nonnegative when sin 9 is nonzero assures 
that sin 0=0 is a saddle point (that is, the direction of a is arbitrary) and not a sin- 
gularity (as contrasted with results of ref. 3). 

In choosing a form for f(0), simplicity is desired to keep the complexity of the 
mechanization at a minimum. The choice given by equation (39), indeed, meets this 
requirement. Now to investigate the behavior of the system for this choice, consider a 
perturbation of the system 

O<|0 Q |«1 ra d (47) 

for no dynamics (p = q = r = 0); then equations (39) to (42) to (44) become 


(48) 

(49) 
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a = -K0 sgn(sin 0) 
0 = -cos (pa 
9 = sin cpa 


(50) 



<P = e \p 


(51) 


Eliminating 0 

and a from equations (49) to (51) gives 



0 = K 0 cos 0 sgn(sin 0 ) 

(52) 


9 = -Kf?|sin 0 1 

(53) 


0 = K 02 C os 0 sgn(sin 0 ) 

(54) 

Note that 

VII 


is small since | 0 | « 1. Therefore assume that 0 (sin <p 0 ± 0) in equations (48), (52), 
and (53) is constant so that the approximate solutions of these equations are, respectively, 

-K|sin 0 o |t 

0 = 0 o e ( 55 ) 


sin 0 O V f 

(56) 


a-a Q + . t ) 

0 sin 0 O x ' 

(57) 


Using equation (55) in equation (54) yields 


* - ^o + -f- 


0n 2 COS <Pr 


sin 0 Q 



It is seen from the form of these approximate solutions that l/^Kjsin assumes 

the role of a relaxation time. The dependence of the relaxation time on sin 0 is not 
a serious problem, provided K is chosen to be sufficiently large, since | sin 0 | is 
"maximized.” 


Since this solution was nondynamical, the space orientation of the two coordinate 
frames is unchanged. The slight perturbation induces motion, in general, in all the gim- 
bal angles such that the orientation is preserved. Hence, the angles are not unique for 
any given orientation. The gain constant K has been related to the relaxation time. 

This relaxation time must be sufficiently small to allow the gimbal servodrives to respond. 
In the following discussion, an upper limit for K is found from this consideration. In 
addition, a lower limit for K is established for the dynamic case to maintain a nonsingu- 
lar solution. 


Let a> L be the maximum obtainable slew rate of the gimbal system. For the non- 
dynamic case, it is seen from equations (48) to (51) that each angular rate is less than or 
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equal to | a |, where 

| ci | = K|sin 0| 

(58) 

Clearly the choice 

«i 

K < L 

(59) 


sin 0 max 



where 0 max is the largest attainable value of | 0| , will limit nondynamically all the 
gimbal angle rates to within the gimbal system tolerances. 

To make a dynamic estimate, consider equation (43). Basic to the derivation of the 
four- gimbal problem is the assumption that a value of a can be found such that cos 0 
always remains large. It is then obvious that the a term must have the dominant mag- 
nitude of all the terms in equation (43) for some 0 not zero; visually 


a sin 0 > |(p sin a + q cos Qf)cos 0 - r sin 0 1 


(60) 


This inequality would be difficult to justify rigorously; however, in view of equation (46), 
inequality (60) is at least plausible. Since the sine and cosine are magnitude bound by 
unity, then 

| p | + |q | + | r | > | (p sin a + q cos a)cos 0 - r sin 0 1 


and 


a = K sin 0 £ a sin 


♦I 


Therefore the choice (for some 0 + 0) 


K|sin 0| > j p | + | q| + | rj (61) 

insures a nonsingular formulation, under the assumption of plausibility of inequality (60). 
Let 

P = sup(|p|,|q|,|r|) 

Then the choice 


K > 


3/3 

si n ^ max 


(62) 


insures inequality (61) for some value of /3 and, hence, insures a nonsingular problem. 
Combining inequalities (59) and (62) yields 

3 £ <K< “L 

| sin 0 m ax | | sin0 max| 

The value 0 max is the upper limit for those values of |0 1 for which inequality (60) 
is not valid; that is, the a term is not sufficiently large to increase cos 0 or, 
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alternately, decrease | 0|. For values of | 9 | on the order of 0 max , the a term is 
dominant and |0| decreases. Hence, 0 max is a limiting value for 1 0 j . 


In choosing a value for 0 max , consideration must be given to the singular point in 
equation (42). The value for 0 max must be sufficiently far from the singular point so 
that excessive rates do not occur. A reasonable choice would be sin 0 max = 0.75 from 
which 


4/3 < K < 


4u>l 

3 


This range for K will maintain a nonsingular solution since 



77 


and the maximum angular rates should not tax the capabilities of the gimbal systems. 


Application 

The problem of computing for 0, 9, 0, and a is formally solved by using equa- 
tion (39), in conjunction with the matrix elements C13, c 21 , c 22 , C33, d23, and d33 
of equations (28) and (33). The resulting equations (with eq. (39) included) are 


a = -K sin 0 sgn(sin 0) 


sin 0 = 


(a 12 sin a + a 22 cos a) 
cos 9 ~ 


(63) 


sin 

cos 

cos 


cos 9 

0 = a 3 j sin 0 - |an sin a + a 2 i cos a ijcos <p 
0 = ( a 12 sin a + a22 cos o:)cos 0 - a32 sin 0 
sin 9 = a 23 sin a - aj3 cos a 
9 = (ai3 sin a + a 23 cos or)sin 0 + a 33 cos 0 


(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 


This solution for 0, 9, 0, and a is algebraic with the exception of a. 

The computation of 0, 9, <p, and a can also be accomplished by solving the 
differential equations (39) and (42) to (44). This solution is equivalent, neglecting com- 
putation error, to the algebraic solution. 
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Mechanization 


Only the mechanization of the algebraic formulation is considered. The matrix 
D(t) has the functional form of the three-gimbal problem previously discussed. Hence, 
the computer circuit in figure 5, with the logic deleted, can be used to solve this matrix. 
From figure 5, the elements 


d 12 = a 12 Ca " a 22 s “ 
d 22 = a 12 Sa + a 22 Ca 
d 32 = a 32 

d 13 = a 13 c “ - a 23 s “ 


(69) 


d 23 = a 13 Sa + ^3°“ 
d 33 " a 33 


are needed to compute the angles \p, 6, and <fi. Since 

sgn(sin <p) = sgn(^) 

a low-accuracy linear or sine potentiometer on the shaft can be used in conjunction 
with dj 3 to compute a. For the linear potentiometer 

a = +Kdj 3 sgn ($) (70) 

Using an integrating servomotor results in the mechanization shown in figure 7. 


a l2 d 12 c 12 S i/i 
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RESULTS AND DISCUSSION 


The computer circuit in figure 5 was programed on an analog computer with A(t) 
generated by the Euler parameter rate equations given in appendix A. 

The Euler angle rate equations admit continuous solutions for constant p, q, and 
r; this is seen from the right- and left-hand limits of equation (13) which has two princi- 
pal solutions at each limit which are equal. Thus, continuous pitch capability is possible 
provided the solution with A 6 = 0 is chosen. The other solution has A<+> = (eq. (14)). 
The logic shown in figure 5 only allows the continuous solutions. If the ratio of q and 
r is changed at the singular points, then discontinuities will result. The logic in figure 5 
minimizes the magnitudes of these discontinuities. 

If a singular point is approached but not encountered, the solutions are continuous 
but the angular rates are exceedingly high. Such a solution is shown in figure 8. This 
solution was obtained by initiating a small roll maneuver (0 = 3°) followed by a pitchover 
(q > 0) and return (q < 0). The solution in this figure was obtained from the circuit in 
figure 5 with the logic disabled. 

If the logic in figure 5 is enabled, then the high rates over the large change in the 
angles t// and 0 , shown in figure 8, are replaced by small discontinuities. The same 
solution is shown in figure 9 with the logic enabled (the dead band was chosen to be 
exceedingly large, ~20°, in order to demonstrate clearly the nature of the solution, 
whereas typically the dead band is chosen to be from 2° to 5°). These high rate solutions 
could result in rate limiting of mechanical equipment which may cause serious degrada- 
tion of the visual display. The small discontinuities are not believed to be as serious a 
limitation because of the relatively small distances over which rate limiting does occur 
(a smaller dead band would decrease the size of the discontinuity shown in fig. 9). 

The method of enhancing 
three-gimbal systems may be 
adequate for many simulation 
problems. When fast gimbal 
systems are available and 
low accuracy is sufficient 
around the singular points, 
a three-gimbal system is 
often suitable, especially in 
view of the cost difference 
between a three- and a four- 
gimbal system. Associated 
with the dead band is a static 



Figure 8.- Solution for pitchover Figure 9.- Approximate solution 

and return with an initial ,or P^ c ^ over an( * return with 

small roll angle. an initial smaM roM an 9 le - 
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error which results in a pointing inaccuracy that is less than or equal to the dead-hand 
size. If a pointing measurement must be made near the singular point, then the dead 
band must be made as small as possible to reduce the static error. However, the 
dynamic error will increase and one must decide from the simulation requirements 
what is acceptable. 

Equations (39) and (63) to (68) for the four-gimbal system were programed on an 
analog computer with the matrix A(t) again generated by the Euler parameter rate 
equations. A pitchover was executed with K = 2, 4, and 6; the results are shown in fig- 
ures 10, 11, and 12, respectively. These solutions have sin 0 = 0 initially, which is 
the saddle point of the inequality constraint. As seen in these figures, the resultant 
motion is entirely absorbed in the 0-axis of the gimbal. Pitchover with sin 0=0 ini- 
tially is where "gimbal flip" should occur according to reference 3. 

It appears from the maximum values of a in figures 10, 11, and 12 that K = 2 is 
adequate when the components of the angular velocity vector are less than or equal to 
1/2 radian per second. Note that this value of K is 4/3, the lower limit established for 
K. With this lower limit, the solution for various inputs p, q, and r is shown in fig- 
ure 13. Again, there is no indication of misbehavior. 

Because Euler angles are needed in many simulation problems, several special- 
purpose analog computers have been constructed with circuits which closely resemble 
that in figure 5. These computers eliminate the task of reprograming the equations for 
the Euler angles of each simulation problem. Furthermore, they conserve the quantity 
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Figure 12.- Pitchover for the four-gimbal model with 
K = 6. 


Figure 13.- Solution of four-gimbal model with K = 2 for 
indicated inputs. 


required of the general purpose simulation computer. Reducing the general purpose 
simulation computer requirements is especially important in digital simulation where 
the resolution of Euler angles often requires as much or more computation time as does 
the remaining portion of the simulation problem. The Euler angle computer is repro- 
grammable for different orders of rotation by simply matching the functional forms of 
the input direction cosines. 

To solve for the four-gimbal angles, the construction of a special-purpose com- 
puter, based on figure 7, is being planned. The computer for four-gimbal systems has 
an added advantage in that it can be used in a three-gimbal mode by bypassing the a 
resolver. Previous all- attitude indicators have been proposed and developed with both 
the three- and four-gimbal concepts used. (See, for example, refs. 3, 4, 9, and 10.) 

For the three-gimbal systems proposed and developed, there is some disagreement 
about the way in which a three-gimbal attitude display should operate. The solutions of 
minimum discontinuity are indicated in references 4 and 9, and those of maximum dis- 
continuity are indicated in reference 10. 

A four-gimbal inertial platform (fig. 14) is proposed in reference 4 so that ’'there 
is no possibility of gimbal lock." The proposed mechanization is to operate with the fol- 
lowing constraint: "A pick-off device, mounted between the first and second gimbal rings 
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at the axis of support of the first gimbal ring, causes the third 
gimbal to be driven in such a manner that the first and second 
gimbal rings maintain an orthogonal relationship during any 
possible aircraft maneuver." This constraint is not to be 
taken literally. It is to be noted that this four-gimbal system 
drives a three-gimbal all-attitude display (8-ball). The pro- 
posed system (according to ref. 4) "lacks sensitivity within 
two or three degrees of the vertical" and thus requires a 
"pitch indicator follow-up system" so that a change in heading 
"is shown as soon as it occurs" in this "exceptional case." 

This same gimbal configuration is again reported in 
reference 3, page 130, where the problems in the exceptional 
case are more clearly discussed: ". . . the fourth gimbal must "flip" or rotate 180 
degrees as pitch angle passes through 90 degrees." 

The four-gimbal systems discussed in references 3 and 4 are analyzed in appen- 
dix C. The difficulty, in all these systems is encountered when 

— - CP -. ff = cos 9 sin a = 0 

30 

which was shown to be a saddle point (and need not be a singular point as indicated by 
ref. 3) if a proper form for & is chosen, namely, 

0 = f(d!)sgn(cos 9 sin a) 

It appears that a slight alteration could greatly enhance these existing four-gimbal 
systems. 



Figure 14.- Four-gimbal system used 
with stable platform. (See ref. 3.) 


CONCLUDING REMARKS 

The technique described in this paper extends the usefulness of the many existing 
three-gimbal systems by providing a method of handling solutions at the singular points 
(gimbal lock) in a predictable and desirable manner. Even so, the quality of motion in 
the neighborhood of the singular points obtainable from a three-gimbal system may not 
be sufficient for all applications. For example, a static error introduced in the solution 
may be as large as the dead band (2° to 5°). 

A four-gimbal system can entirely eliminate discontinuities. The experimental 
results from the four-gimbal model indicate no adverse behavior. Due to the non- 
uniqueness of the angles, error in the computations of the fourth angle does not affect the 
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static or dynamic accuracy of the gimbal system when the essentially algebraic formula- 
tion is used. 

It appears that an improved version of a four- gimbal stabilized platform could be 
built by applying the results of this investigation. Clearly, the need for a feasibility study 
is indicated by the analysis. 

The two special-purpose analog computers for solving the three- and four-gimbal 
problems provide a cost advantage by reducing programing time and the demands placed 
on general-purpose computing equipment. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 14, 1968, 

125-19-06-02-23. 
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APPENDIX A 


THREE -GIMBAL MATHEMATICAL DETAIL 


The matrix A(t) is found- (ref. 8) as a function of the Euler parameters (a normal- 
ized quaternion) by 


where 


A(t) = 


2(a2 + d2) - 1 2(ab - cd) 2(ac + bd) 

2(ab + cd) 2(b2 + d2) - 1 2(cb - ad) 

2(ac - bd) 2(cb + ad) 2(c2 + d2) - 1 


a = i(-pd - qc + rb) 


b = |(pc - qd - ra) 
c = |(-pb + qa - rd) 
d = ^-(pa + qb + rc) 


(Al) 


(A2) 


and p, q, and r are the components of the angular velocity vector measured along the 
rotating coordinate axis. 


A <t) = M x [«] y M z = 


of the three Euler angles 

xf/, 9, and 

(c\pc e) 

(S\pC6) 

(-S B) ~ 

f tipsesA 

(CilsC<t> \ 

(s<pce) 

\-s*c<p ) 

\+S'(pSxl/SdJ 


/S^s <t> \ 

fC(pS\pse\ 

(c<pce) 

\+c<pC\is sej 

[-S^Cxp ) 



(A3) 
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APPENDIX A 


Provided $ is known, a matrix B(t) can be computed as 


^11 a 12 

B (t) = = ( a 21 C< ^ ~ a 31 s <£>) ( a 22 c <t> ~ a 32 s <f) 

( a 21 s 0 + a 31 C0) (a22 s 0 + a 32 c 0) 


( a 23 C0 - a 33 S0) 
( a 2 3 S^ + a 33 C^ 


(A4) 


where B(t) is functionally related to 9 and t// by 

~C\l/Ce S1//C0 -so, 

B (f) = 00yK1 z = " S ^ 0 (A5) 

Ci//S0 s^se ce 
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APPENDIX B 


FOUR- GIMBAL MATHEMATICAL DETAIL 


Expanding A(t) as a function of the four gimbal angles \p, 0, cp, and a yields 

/CaCxjsCd \ /CaS^Ce \ f-S9Ca \ 

[+SaCi//S0S4>] l+SaC\pC<p \+SaC0S<p) 

\-SaS\pC<t> ) \-SaS(l>Sil/S^ 


A(t) = [a] £ Mx ^ Y ^ Z 


/CotCxpSdSA /CaCxjy Cep \ /S0Sa \ 

J -CaStl/Ccp ( +CaS<J>S<//S0 \+S0C0 Cqi/ 

\-SaC\f/C9 ) \-SaS\jyCe J 


/SipScp \ (C<pSipS9\ 

\+C(pC\i/Sd) \-s^Ci// / 


With a knowledge of a, the matrix D(t) is 


D(t) = [-q] z A(t) = 


( a llCQ! - a 2 iSa) 
(auSa + a^Ca) 


l 31 


(a 12 Ca - a 22 Sa) 
(a 12 Sa + a 22 Co!) 

a 32 


(C0C0) 


( a 13 C« ' a 23 Sc*) 
( a 13 S« + a 23 Co) 
a 33 


The matrix D(t) as a function of the three angles i//, 0, and <p is given by 


(Bl) 


(B2) 


D(t) = [0 ] x [«] y M z = 


(c^C0) (S\]/ce) 

/CI//S0S <A /Cil/Ccp \ 

[-Sxpccp ) \+S4>Sxj/Sel 

fsxj/S<f) \ fc<ps\pse\ 

\+C(j)C\l/Sd) \-S 0 C^ J 


(-S0) 

(s<pce) 


(C(pce) 


The solution for \p, 0, and 4> from equation (B3) is well known. 


(B3) 
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APPENDIX B 


Again, C(t) is found to be a function of A(t), a, and <p as 


C(t) = AW = 


(a n Ca - a 21 Sa) 

/a^C^S a + a.2\C<t)C(A 
\" a 31 S< ^ / 

/aiiS^Sa + a 2 ^S0Co' 
\+ a 31 c 0 


(a 12 Ca - a 22 Sor) 

/a^Cc^Sof + a 22^ < ^^' a \ 
\ _a 32 s< ^ / 

/a^S^Sa + a^S^Ca 1 

l +a 32 C< ^ 


( a l 3 Ca - a 23 Sa) 

/a^Ct^Sa + a23C0CcA 
\- a 33 s 4> / 

(a^S^Sa + a22®0C^ 
+a 33 C^ 


(B4) 


or of 6 and as 



Ci//Ce 

s\p>C6 

-S0 

c(t> = m y bo z = 

-Si^ 


0 


C1//S0 

S0/S0 

C0 


(B5) 


29 



appendix c 


ANALYSIS OF SOME EXISTING FLIGHT SYSTEMS 


For the four-gimbal configuration shown in figure 14, the transformation from the 
fixed to the rotating frame is given by 

V r = M x [o] Y [«] x M z Vf = A(t)Vf (Cl) 

The angle cp must be chosen such that a is not allowed to approach ±90° (that is, 
cos a must not be zero). Hence, the necessary constraint is 

8cos q & > o (c2) 

d<p 

where $ must not be zero whenever cos a differs from 1. The expression of cos a 
as a function of A(t), <//, 9, and <p is found from the following matrix equation: 

B (t) = [«] X M Z = [- 0 ] Y Etl x A(t) (C3) 

from which 

b 3 g = cos a = -a 13 sin 9 + cos 9(sl 23 sin cf> + a 33 cos ^ (C4) 

The partial derivative in equation (C2) is then 


Use is then made of 


9 cos a 


= cos 0^23 cos (p - a33 sin <p^ 


b 3 3 = sin a = a 3 3 cos <p - a 33 sin cp 


from which 


9 cos a 

dtp 


= cos 9 sin a 


The choice for cp is then 

cp = f(o!)sgn(cos 6 sin a) 

where i(a) is a positive bounded function that is nonzero for cos a * 1. 
The form for cp as indicated by references 3 and 4 is 

4> = f(a) sgn(o?) 


(C5) 


(C6) 


(C7) 


(C8) 


which is insufficient to maintain a nonsingular formulation. Note that "gimbal flip” and 
the "exceptional case" occur when cos 9 changes sign. 
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